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Abstract 

We apply the formulation developed in a recent paper [Y. Ohashi and A. Griffin, Phys. Rev. 
A 72, 013601, (2005)] for single-particle excitations in the BCS-BEC crossover to the case of a 
broad Feshbach resonance. At T = 0, we solve the Bogoliubov-de Gennes coupled equations taking 
into account a Bose condensate of bound states (molecules). In the case of a broad resonance, 
the density profile n(r), as well as the profile of the superfluid order parameter A(r), are spatially 
spread out to the Thomas-Fermi radius, even in the crossover region. This order parameter A(r) 
suppresses the effects of low-energy Andreev bound states on the rf-tunneling current. As a result, 
the peak energy in the rf-spectrum is found to occur at an energy equal to the superfluid order 
parameter A(r = 0) at the center of the trap, in contrast to the case of a narrow resonance, and 
in agreement with recent measurements. The LDA is found to give a good approximation for the 
rf-tunneling spectrum. 

PACS numbers: 03.75.Ss, 03.75.Kk, 03.70.+k 
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I. INTRODUCTION 

In a recent paper we presented a microscopic theory (at T = 0) of the single-particle 
excitations in the BCS-BEC crossover regime of a trapped superfluid Fermi gas with a Fes- 
hbach resonance. This theory included a Feshbach resonance and the associated molecules. 
The Bogoliubov single-particle excitations were found by numerically solving the microscopic 
Bogoliubov-de Gennes (BdG) coupled equations^]. The calculations also involved calculat- 
ing self-consistently the equation of state given by the Fermi chemical potential (which takes 
into account the depletion of the number of Fermi atoms as they form into stable molecules). 
In Ref. jll, we applied this theory to a narrow Feshbach resonance {gr^/n < e-p, where (?r 
is the coupling constant of the Feshbach resonance, n is the number density of atoms, and 
Ep is the Fermi energy of a non- interacting Fermi gas). We calculated in a self-consistent 
manner the density profile, superfluid order parameter, and the single-particle excitation 
gap Eg in the entire BCS-BEC crossover region. We also discussed the effect of Andreev 
bound states[3] on the rf-tunneling current dataj^. 

In this paper, we extend our previous work and present results for single-particle Fermi 
excitations in the BCS-BEC crossover regime in the case of a broad Feshbach resonance 
[gry/n ^ e-p]. This extension is veryrelevant because all current experiments on superfluid 
Fermi gases ^°K|5:] and ^Lij^, Q, 0,0] make use of a broad Feshbach resonance in order to 
tune the pairing interaction. In the case of a broad Feshbach resonance, the BCS-BEC 
crossover occurs when the bare threshold energy 2z/ of the Feshbach resonance is much 
larger than the Fermi energy ep, and as a result, the Fermi atoms in the open channel 
dominate. The Feshbach resonance contributes to an effective pairing interaction between 
Fermi atoms in the open channel, but the number of stable Feshbach molecules is very small 
in the crossover region for a broad Feshbach resonance. This situation is quite different from 
a narrow resonance case, where the Feshbach molecules already dominate in the crossover 
region. In this paper, we will show how single-particle BCS-like excitations given by the 
BdG equations differ for the two types of resonances. In particular, we will discuss how this 
difference shows up in the low-energy spectrum of rf-tunneling current dataj^. 

We stress that while the differences between the results we find for the excitation spectrum 
in the case of a broad and narrow resonance are substantial, they emerge from the same 
microscopic theory based on the Bogoliubov-de Gennes coupled equations. Ultimately the 
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differences arise from quantitative changes in the BdG solutions, such as the spatial width 
of the order parameter in a trap. 

The numerical solutions of the BdG equations for a broad resonance which we report 
here are computationally very challenging, because of the need to expand the basis of states 
to much higher energy. We briefly discuss the associated cutoff problems in Sec. II. 

As has been discussed extensively in the recent literature, a broad Feshbach resonance can 
be treated using the single-channel model, pioneered by Leggett jl^. Nozieres and Schmitt- 



Rink[ll|, RanderiaQ, 

and coworkers |l3l. Ref.|l| and the present paper are built on 

a coupled fermion-boson two-channel model (see Sec. IT). This is necessary to deal with a 
narrow Feshbach resonance, which was treated in Ref.|l|. We think that treating a broad 
Feshbach resonance using a two-channel model, while not necessary, has an advantage in 
that it allows one to see how the difference between a broad and narrow Feshbach resonance 
emerge within a unified formalism. However, we note that within a single- channel model, 
our BdG equation approach can be viewed as an extension of the original work by Bruunjlfij]. 
which only considered the BCS limit of the BCS-BEC crossover region. 

This paper is organized as follows. Since the formalism presented in Ref. Ilj] is valid for 
both narrow and broad Feshbach resonances, we refer to this earlier paper for most details. In 
Sec. II, we briefly summarize the coupled fermion-boson model, which our analysis is based 
upon. The proper treatment of higher energy states in the numerical calculations for the case 
of a broad Feshbach resonance is noted. In Sec. Ill, we present our self-consistent numerical 
solutions of the BdG equations for a broad Feshbach resonance. We discuss single-particle 
excitations in Sec. IV and use these results to calculate the rf-tunneling current response in 
Sec. V. 



II. COUPLED FERMION-BOSON MODEL 



A superfluid Fermi gas with a Feshbach resonance can be described by the coupled 
fermion-boson model 

H = J2j t^rvl/t(r)[-^-/i + V;f,p(r)]v^.(r)-f/| «{(r)vl/t(r)vl;^(r)v^Kr) 



+ / (ir<l>^(r) 



2M 



+ 2z/-/XM + Cap(r)l*(r)+^7. / rfr [<|.t(r)vl/^(r)vl/^(r) + /i.c 



(2.1) 
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We refer to Ref. [l| for further discussion about this model. We assume a two-component 
Fermi gas described by a fermion field operator ^l^o-(r), where the pseudo-spin label cr =1, | 
represent two atomic hyperfine states. These states are in the so-called open-channel. $(r) 
describes the molecular bosons associated with a Feshbach resonance, is a coupling 
constant describing the Feshbach resonance, and 2v is the bare threshold energy of the 
resonance. Experimentally, 2v is tunable by an external magnetic field, so that the effective 
pairing associated with the Feshbach resonance is also tunable. This interaction becomes 
stronger as we decrease the threshold energy 2v. The atoms and molecules are trapped in 
harmonic potentials V^^^ir) = {1 j 2)mujQr'^ and V^^-p{r) = {l/2)Mu!Qr'^, respectively. Here, 
we assume an isotropic trap, for simplicity, and assume that atoms and molecules feel the 
same trap frequency c^o (as is the case in current experiments). 

Since a molecule consists of two Fermi atoms, the molecular mass is twice as large as a 
atomic mass, M = 2m. We also impose the conservation of the total number of atoms, 

N = Nf + 2Nm = J rfmF(r) + 2 J dmM(r), (2.2) 

where A^f and A^m represent the number of Fermi atoms in the open channel and the number 
of molecules associated with the Feshbach resonance, respectively. Equation ()2.2j) defines 
the number density of Fermi atoms in the open channel np{r) and the local density ny[{r) 
of molecules. Since we limit ourselves to zero temperature in this paper, all the molecules 
are Bose-condensed. The constraint in Eq. ()2.2j) has already been taken into account in 
Eq. ()2.H) by taking the Bose chemical potential fiM = '^fJ', where fi is the Fermi chemical 
potential. Equation ()2.1|) also includes a non-resonance interaction U, which is taken to be 
weakly attractive. 

The BCS-BEC crossover at T = is treated in Ref. by solving the BdG coupled equa- 
tions for the mean-field approximation to the Hamiltonian in Eq. ()2.H) . together with the 
equation for the number of atoms |lfl l|. U sing Eq . (| 2.1D. the superfluid phase is characterized 
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by the composite order parameter! 

A(r) = A(r)-(7,0Af(r). (2.3) 

Here, A(r) = f/(\&|(r)^'|(r)) is the Cooper-pair order parameter, and 0a/ = (^(r)) describes 
ecular condensate. For details about the self-consistent calculation of A(r), we refer 



the mo 
to Ref. 
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In this paper, we take g,. = g^yN/Rp = lOep to describe a broad Feshbach resonance, 
where Rp = 2eF / mujQ is the Thomas-Fermi radius. This is the only parameter change from 
Ref. Q], which considered a narrow resonance with g^ = 0.2eF. The total number of atoms 
is again taken to be = 10912, which corresponds to the Fermi energy ep = SLScuq, where 
ujq is the trap frequency. For the non-resonant interaction U, we take [/{N/R^) = O.SSep- 

In solving the BdG equations numerically, one cannot retain all the harmonic potential 
eigenfunctions. Thus we need to truncate the basis states by introducing a finite cutoff 
energy lUc- In this paper, we take Uc = IGLScJo = 5.1£:f- To avoid any spurious effects of this 
cutoff, we need to restrict our calculations to the region where the maximum value of the 
composite order parameter A(r) in the center of the trap is much smaller than this cutoff 
Uc- In the case of a narrow Feshbach resonance {g^ ^Jn/R^ < ep), A(r) does not become 
very large even in the BEG regime, so that this computational restriction is not a serious 
problem. However, A(r) in a broad Feshbach resonance (g-^ \J N/ Rp ^ ef) can be very large 
in the BEG regime. Thus, in this paper, we restrict our calculations in the strong-coupling 
BEG regime to the region (/cpag)^^ < 1, where A(r = 0) ~ 0(6^) (here kp is the Fermi 
momentum, and is the two-body atomic s-wave scattering length). In the equation of 
state involving the number of atoms, we still take into account the contribution coming 
from the high energy states with energies u > hj using a local density approximation 
(LDA). These high energy contributions are important, since neglecting them means that 
the number of molecules um is overestimated in the BEG regime. 



III. DENSITY PROFILE AND COMPOSITE ORDER PARAMETER 

Figure 1 shows our results for the Fermi chemical potential /i in the BGS-BEG crossover. 
This is very similar to the result for a narrow Feshbach resonance (see Fig. 6 in Ref. 
^1). That is, for both broad and narrow Feshbach resonances, the chemical potential 
has qualitatively the same dependence on the value of the parameter (/^pfls)^^, where as 
is the two-body scattering length. We note that fi gradually deviates from being equal to 
the Fermi energy ep = 31.5c(Jo as we increase the pairing interaction, and changes sign at 
(/cpcts)^^ — 0.65. For a broad Feshbach resonance, this change is found to correspond to 
u ~ SOep (see the inset of Fig. 1). This threshold energy is much higher than in the case 
of a narrow resonance (see Fig. 6 of Ref.|l|), where we found that yU = at 2z/ ~ 0. For 
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a broad resonance, we might note that the generahzed s-wave scattering length a<j defined 
in Eq. (4.5) of Ref.jll reduces to the usual two-body scattering length of for a Feshbach 
resonance defined in Eq. (4.7) of Ref.jll. This is because we have 2v ^ 2/i in the crossover 
region for a broad resonance. 

Figure 2 shows the total atomic density profile n{r) = nF{r) + 2nM{r). In the BCS 
regime, the density profile is significant right up to the Thomas- Fermi radium Rp, which 
describes the size of a non-interacting Fermi gas in a harmonic trap at T = 0. The density 
profile gradually shrinks as one approaches the (strong-coupling) BEC regime. However, 
this narrowing is much less than we found in the case of a narrow resonance (see Fig. 5 in 
Ref. j.l|). In a broad resonance, the Fermi atoms in the open channel are still the dominant 
particles in the crossover region. (As noted above, this occurs at large values of 2z/, as shown 
by the inset in Fig. 1.) As a result, repulsion resulting from the Pauli exclusion principle 
between Fermi atoms gives a wider atomic density profile in the case of a broad Feshbach 
resonance. This is quite different in the case of a narrow resonance, where the Feshbach 
molecules (bosons) already dominate in the crossover region, and the Pauli principle is less 
effective. As a result, the spatial profile of n(r) is more localized at the center of the trap 
in a narrow resonance. 

Although the Fermi atoms in the open channel are dominant in the broad Feshbach 
resonance, we note that even at T = 0, only a small fraction of the atoms are associated with 
Cooper-pairs. Figure 3 shows the Cooper-pair condensate fraction A^,^^, defined by[2fl[ I21I1 



In the BCS regime. Fig. 3 shows that only a small fraction of A^f contributes to the Cooper- 
pair condensate fraction A^^ . This refiects the fact that only atoms near the Fermi surface 
form a Cooper-pair Bose-condensate in the weak-coupling BCS regime. Atomic states well 
below the Fermi energy are essentially unmodified from those in a non-interacting Fermi 
gas. The Cooper-pair condensate fraction starts to increase as one approaches the BEC 
regime. At (kpas)^^ ~ 1, we find 2Nq ~ Np, indicating that most atoms in the open 
channel form Cooper-pairs which are Bose-condensed. Namely, this regime can be regarded 
as the BEC of bound pairs of atoms from the open channel. 

Figure 4 shows the spatial variation of the composite order parameter A(r) in the 
crossover region. Although the number of Feshbach molecules A^m is very small (see Fig. 3), 




(3.1) 
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we find that the molecular BEC order parameter —gr<pM{f^) is comparable to the Cooper- 
pair component A(r). Since the coupling is large in a broad Feshbach resonance, the 
contribution of the molecular order parameter —gr4>M{f') is clearly enhanced. Even in the 
BCS region at (fepa..)"^ = —1 [where Nm/N^ ~ 10^^, as in Fig.3], Fig 4(a) shows that 
— (7i.$m(0)/A(0) ~ 0.24. In Fig. 4(d), the molecular BEC component becomes dominant 
over the Cooper-pair component, although the number of condensed Feshbach molecules is 
still very small, as shown in Figs. 2(d) and 3. Thus, even in the case of a broad Feshbach reso- 
nance, the Feshbach molecules are found to make a major contribution to the composite order 
parameter A(r), although the number of the Feshbach molecules is very small. This surpris- 
ing result is due to the fact that the Cooper-pair order parameter A(r) = (\l/|(r)\I'|(r)) is 
quite different from the number density of atoms np{r) = J^ai'^ti''^)'^ i'^ ^ Fermi super- 
fluid. The situation is quite different from a Bose gas BEC, where there is a simple relation 
between the BEC order parameter 0Af('") = {^{f^)) and the number density of condensed 
bosons at T = 0, namely UM^r) = ($"l'(r)$(r)) = |0M(r)p. The fact that nM_{r)/np{r) <^ 1 
does not necessarily mean a small value of the ratio — (7r0M('")/A(r). 



IV. SINGLE-PARTICLE EXCITATION SPECTRUM 

Figure 5 shows the fermion single-particle local density of states, defined byQ] 

N{uj, r) = --Im[G'ii(r, r, ^ + i5)]. (4.1) 

Here Gn is the analytic-continued BCS-like single-particle thermal Green's function, 

G'n(r,r',z^„) = - dre'"^-^ {Tr{^^{r,T)^{r' ,0)}) , (4.2) 
Jo ' 

where Un is the fermion Matsubara frequency. Equation 1)4.111 describes single-particle exci- 
tations of energy uj at position r. 

As discussed in Ref. atoms in a trapped Fermi superfluid feel the effect of a combined 
potential well, consisting of a diagonal trap potential V^f^^^r) — fi and an off-diagonal pair 
potential A(r). In the BCS regime, this combined potential well has a minimum at r ~ Rp 

nri 

(where Rp is the Thomas- Fermi radius). The Andreev bound states [1|, y] are then formed 
near the bottom of this potential well, which is near the edge of the trapped gas (~ -Rf)- 
Indeed, we find some Andreev bound states localized around r ~ i?F in Fig. 5(a), and the 
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lowest one gives the single-particle energy gap Eg (~ uq) of this trapped Fermi superfluid. We 
note that Eg is much smaller than the typical magnitude of the composite order parameter 
at the center of the trap. [In Fig. 5(a), we find that A(r = 0) ~ IOujq ^ Eg.]. In Fig. 
5(a), we see that there is a large excitation gap (~ IOujq ^ Eg) near the center of the 
trap (r ~ 0). This large energy gap in N{uj,r) is seen to be equal to the magnitude of the 
composite order parameter A(r = 0) in the center of the trap. Thus while the true threshold 
energy of the single-particle excitation spectrum is determined by the lowest Andreev states 
at the edge of the trap, the local excitation gap at the center of the trap is equal to the 
value of the local composite order parameter A(r = 0). In the BCS regime, since Eg is 
largely determined by the width of the combined potential well. Eg is not very sensitive to 
the increase of the composite order parameter A(r = 0) in the center of the trap. As a 
result, in the BCS regime. Eg only slowly increases as we go from the BCS region to the 
crossover regime (where {kpas)~^ ~ 0). This is to be compared with the rapid increase of 
A(r = 0) on the BEC side, as shown in Fig. 6. 

In the BEC regime, when the Fermi chemical potential /i is negative. Fig. 5(b) shows 
that the local density of states N{uj,r) has a large energy gap everywhere, including at 
the edge of the trap. The physics of this can be easily understood, using the Bogoliubov 
single-particle excitation spectrum for a uniform Fermi superfluid given by 



When /i < 0, this energy spectrum Ep has an energy gap|10] given hy Eg = \j jj? + A"^. In 
the (strong-coupling) BEC limit, this gap approaches Namely, instead of the superfluid 
order parameter A, it is the magnitude of the Fermi chemical potential which determines 
the single-particle excitation gap. In turn, this determines the binding energy of a Cooper 
pair. Similarly, in a trapped gas, even though the order parameter is small near the edge 
of the gas. Eg is still large in the BEC regime due to a large magnitude of the chemical 
potential Figure 6 shows how the calculated energy gap Eg approaches in the BEC 
regime (/x < 0), as expected. ^ 

In a narrow Feshbach resonance considered in Ref. p|, the single-particle excitation gap 
Eg has a minimum value when one is in the crossover region (see Fig. 11 of Ref. 0]). This 
is due to the significant narrowing of the order parameter spatial profile in the crossover 
region. For a narrow resonance, the width of the combined potential well becomes wider 




(4.3) 
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in the crossover region, and consequently the energy of the lowest Andreev state decreases. 
On the other hand, for a broad Feshbach resonance, the narrowing of the order parameter 
spatial profile is not significant (see Fig. 4). That is, the width of the combined potential well 
does not become wider, in contrast with the case of a narrow resonance. As a result, the 
single-particle energy gap Eg, determined by the lowest Andreev state, steadily increases 
as we move from the BCS to BEG region. However, this increase of Eg is still slow in 
the BCS regime, compared with the rapid increase in A(r = 0), as shown in Fig. 6. In 
the next section, we will point out that this quite different behavior of Eg in the case of 
a broad and narrow Feshbach resonance has important implications for the rf-tunneling 
current spectroscopy. 

V. RF-TUNNELING CURRENT SPECTRUM 

In Figs. 7 and 8, the solid line shows the calculated rf-tunneling current spectrum If{uj) 
in the BCS-BEC crossover region. Here u = ujl — uja — + l^a is the effective detuning 
frequency. The frequency of the laser light is ujl-, while Ua and /i^ are the threshold energy 
and chemical potential, respectively, of the additional atomic hyperfine state \a) which the 
open state | is coupled to. For the detailed derivation and discussion of the expression for 
If{uj) we use, based on the solution of the microscopic BdG equations, see Sec. VII of Ref. 

Q. 

In a broad Feshbach resonance, we find that the peak energy in the rf-spectrum increases 
with increasing pairing interaction [panels (a)-(c)]. In the BEC regime [panel (d)], the 
spectrum /f(co') also exhibits a finite energy gap. This energy gap reflects the large single- 
particle excitation gap in the BEC regime, shown in Fig. 5(b) and discussed in Sec. IV. 

In Fig. 6, the solid circles show the peak energies of the calculated rf-tunneling current 
spectra in the crossover region (as shown by the examples in Figs. 7 and 8). These T = 
results clearly indicate that for a broad resonance, the peak energy in the rf-spectrum can 
give a direct measurement of the magnitude of the order parameter A(r = 0) in the center of 
the trap. Figure 6 also shows the ohserved'peak energies in the recent rf-tunneling experiment 
on superfluid ^Li, as open circles j4^. Although the present calculation is for an isotropic trap 
at T = while the experiments have been done in a cigar-shaped trap at finite temperatures, 
we note that our results agree reasonably well with the experimental data, and especially 
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well in the BEC region [i.e., in the region {kpas)~^ > 0.5]. 

In the case of a narrow Feshbach resonance the very low-energy rf-spectrum in the 
BCS regime is dominated by Andreev excitations localized near the edge of the trap. As 
discussed in detail in Sec. VII of Ref.|l!], the result is that the contribution of the large 
order parameter in the center of the trap is hidden in If{uj) by the spectral weight of the 
low-energy surface excitations. This effect is still significant in the crossover region. This is 
because the spatial profile of the composite order parameter A(r) narrows in the crossover 
region, as shown in Fig. 5 of Ref . Q . As we have discussed in Sec. IV above, this narrowing 
lowers the energies of the Andreev bound states localized at the bottom of the combined 
potential well. Thus for a narrow resonance, these low-energy Andreev states dominate 
the quasi-particle excitation spectrum, with the result that the low-frequency spectrum of 
If{uj) is also dominated by these surface excitations. In the BEC regime (where fi < 0), 
we recall that the effects of these low-energy Andreev-like excitations are suppressed. In a 
broad Feshbach resonance, in which the spatial narrowing of the profile of A(r) was found 
to be not significant, the role of the low-energy quasi-particle excitations at the edge of the 
trap is less dominant, compared with the narrow resonance case. Our results thus verify [2^ 
that for a broad Feshbach resonance, the peak in the rf-tunneling current spectrum gives a 
direct measurement of the order parameter at the center of the trap. The effect of the low 
energy Andreev bound states is negligible on If{uj). 

In principle, one could try to observe the true single-particle energy gap Eg from the 
threshold energy of the rf-tunneling spectrum. For example, in Fig. 7(b), Eg is given by 
the lowest peak in the plotted rf-spectrum. Since Eg is very small, high resolution would 
be necessary to observe this small gap structure in the rf-tunneling spectrum. However, we 
stress that the these Andreev surface states are of great interest, representing as they do a 
very interesting aspect of the quasiparticle spectrum of a Fermi superfiuid. 

In Figs. 7 and 8, the dashed line shows our results using for If^uj) calculated in the 
local density approximation (LDA). (The LDA expression for If{uj) is given by Eq. (7.19) 
in Ref.y|.) If{^^) was first computed using a LDA in Ref. Q]. Our LDA calculation is 
somewhat different, since self-consistent BdG solutions for A(r) and fi are used. We find 
that when the self-consistent solution of the BdG equations are used, the LDA (dashed 
lines) gives a good overall approximation to the fully microscopic calculations (solid lines) 
in the crossover region. In particular, the peak energy in the LDA result for /f(^^), which 
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corresponds to the magnitude of the composite order parameter A(r = 0) in the center of 
the trap, is in excellent agreement with the microscopic result based on the BdG eigenstates. 



VI. SUMMARY 

In this paper, we have extended our previous theoretical calculations for a narrow Fes- 
hbach resonance to the case of a broad Feshbach resonance. This extension is important 
because all current experiments make use of a broad Feshbach resonance. We have numeri- 
cally solved the BdG coupled equations, together with the equation of state for the number 
of atoms, in the BCS-BEC crossover region at T = 0. This formalism we use was developed 
in detail in Ref. Q]. The case of a broad resonance involves much more demanding numerical 
calculations because of the need to include high energy states in solving the BdG equations. 



For a broad Feshbach resonance {g^-y N/ Rp ^ ep), the BCS-BEC crossover occurs in the 
region where the threshold energy 2z/ of the Feshbach resonance is much larger than the 
Fermi energy ep (see Fig. 1). As a result, molecules associated with the Feshbach resonance 
only contribute to an effective paring interaction between atoms in the open channel in a 
virtual sense. As expected, the number of Feshbach molecules is found to be very small 
compared to the number of Cooper-pairs. Because of this, the effect of Pauli exclusion 
principle between Fermi atoms is still important in the crossover region. As a result, the 
spatial profile of the composite order parameter A(r) is spread out and broad even in the 
crossover region. This is quite different from the narrow resonance case ^ , where Feshbach 
molecules dominate in the crossover region. In this case, A(r) is spatially more centered in 
the trap. 

This difference in the width of the spatial profile of the composite order parameter A(r) 
between a broad and narrow Feshbach resonances leads to qualitatively different results 
for the low-energy rf-tunneling current spectrum. For a broad Feshbach resonance, the 
narrowing of the spatial profile of A(r) is not significant in the crossover region. As a result, 
the energies of the Andreev bound states localized at the edge of the trap are not lowered 
and these states end up having little weight in the rf-tunneling current spectrum If{u!). 
As a result, the peak energy in the rf-spectrum is predicted to occur at the value of the 
composite order parameter A(r = 0) at the center of the trap. In contrast, in a narrow 
Feshbach resonance, the Andreev bound states strongly modify the low-energy rf-tunneling 
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current spectrum in the crossover region. This hides the appearance of the peak at A(r = 0) 
in Ipi^^), as discussed in Ref. 

In summary, for a broad resonance, we have solved the BdG equations for the coupled 
fermion-boson model in Eq. ()2.H) at T = 0. We have taken into account the formation of 
a composite order parameter composed of both Cooper-pairs and a BEC of real molecules. 
The formalism we use was developed in Ref . , where it was applied to the case of a narrow 
Feshbach resonance. In this paper, we used our BdG eigenstates to calculate the rf-tunneling 
current and found excellent overall agreement with the earlier work based on a simpler 
LDaIs^. In particular, we confirm at a more microscopic level (at T = 0) that the finite 
energy peak[4] which develops in the rf-tunneling current in the superfluid phase can be 
used to give a direct measurement of the magnitude of the order parameter at the center of 
the trap. 
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1 0.8 0.6 0.4 0.2 -0.2 -0.4 

(kpas)' 

FIG. 1: Fermi chemical potential // in the BCS-BEC crossover region for a broad Feshbach 
™ce fev^. 10..)^ T,.e i„.„c«on i. .ive„ i„ ...s of .he .w.bod. .wave atomic 
scattering length as\^. The inset shows the same chemical potential // plotted as a function of the 
bare threshold energy 2iy of the Feshbach resonance. 
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FIG. 2: Total atomic density profile n(r) in the BCS-BEC crossover for a broad resonance. We 
also show the atomic density profile in the open channel nF{r) = J2a{^ai^)^'^i^))^ ^ ^^^^ 
the molecular condensate density associated with the Feshbach resonance nuir) = |(/>m(i')P- The 
molecular condensate is very small in panels (a)-(c) [nM(r) ^ n{r)\. The inset in panel (b) shows 
the molecular density profile in more detail. 
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FIG. 3: Character of particles at T = in the crossover region in the case of a broad Feshbach 
resonance. iVp is the total number of Fermi atoms in the open channel, and A^m is the molecular 
condensate fraction (i.e., the number of Bose-condensed Feshbach molecules). Nq is the Cooper- 
pair condensate fraction, describing the BCS phase. 
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FIG. 4: Profile of tlie composite order parameter A(r) in the BCS-BEC crossover for a broad 
resonance. Tlie separate contributions from the Cooper-pair component A(r) and the molecular 
BEC order parameter component —gr(pM{r) are also shown. Compare with the analogous results 
for a narrow resonance shown in Fig. 7 of Ref.[]|. 



18 



N(a),r) [cOo^Rf] v/£f = 250, (kpas)"^ = -0.48 




FIG. 5: Local single-particle density of states N{u;,r). (a) BCS region, (b) BEG region. We 
have introduced a small imaginary part F = 0.2a;o to the energies given by the BdG equations, to 
smooth out the results. 



19 




FIG. 6: Single-particle excitation gap Eg in the BCS-like density of states N(uj,r), and the 
composite order parameter A(r = 0) at the center of the trap. Eg becomes large in the BEG 
region where the chemical potential /i is negative, and it approaches (dotted line) in this strong- 
coupling regime. The solid circles show the highest peak energies in the calculated rf-tunneling 
current spectra Ipi^) shown in Figs. 7 and 8. The open circles show the measured peak energy 

" n 



in the rf-tunneling spectrum in superfluid Li 



for comparison. See Fig. 21 of Ref._l] for the 



analogous results for a narrow Feshbach resonance. 
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(a) v/£f = 600, (kpas) '= - 1 




(b) v/ep = 200, (kpa,) '= - 0.32 
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CO/Ep 

FIG. 7: Calculated rf-tunneling current spectrum in the BCS and crossover regions, based on the 
self-consistent solution of the BdG equations. The solid line shows the result using the microscopic 
expression for the rf-tunneling current spectrum, given by Eq. (7.10) of Ref.Q|. The fine structure 
in the spectrum originates from the discrete BdG excitation energies in the harmonic trap. The 

n 

dashed line is based on the LDA expression, given by Eq. (7.18) of Ref. 1]. For clarity, we have 
introduced a small imaginary part F = 0.5a;o (damping) to the BdG eigenstates to smooth the 
results. 
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FIG. 8: Same plot as in Fig. 7, for (fepas) ^ = 0.66 (the BEC regime). 
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